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ABSTRACT 


Theoretical  investigations  for  fully  developed  laminar 
heat  transfer  in  noncircular  ducts  were  carried  out  using 
both  the  point  matching  and  variational  methods.  The  ther¬ 
mal  boundary  conditions  of  axially  uniform  wall  heat  flux 
and  equal  wall  temperature  at  a  given  axial  position  are 
considered  for:  (a)  circular  ducts  with  diametrically  op¬ 
posite  flat  sides:  (b)  ducts  of  multiply  connected  cross 
sections  in  the  form  of  regular  polygons  with  a  central 
circular  hole  and  circular  ducts  with  a  central  regular 
polygonal  hole.  Heat  transfer  results  for  regular  poly¬ 
gonal  ducts  with  thermal  boundary  conditions  of  both  peri¬ 
pherally  and  axially  uniform  heat  flux  were  obtained  using 
the  point  matching  method. 

Extensive  numerical  results  are  presented  for  local 
quantities  such  as  the  velocity  distribution,  temperature 
distribution,  temperature  gradient  at  the  wall  and  for 
overall  quantities  such  as  the  mass  flow  rate,  friction 
factor,  average  wall  temperature,  bulk  fluid  temperature, 
and  Nusselt  number.  Graphical  results  are  presented  for 
the  shear  stress  distribution,  temperature  gradient  at 
the  wall,  mass  flow  rate,  friction  factor  -  Reynolds 
number  results  and  average  Nusselt  number.  Velocity  and 
temperature  profiles  are  also  presented  for  the  represent¬ 


ative  cases. 
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CHAPTER  I 


INTRODUCTION 


Laminar  forced  convection  heat  transfer  in  non¬ 
circular  ducts,  under  developing  or  fully  developed 
conditions,  has  been  the  subject  of  investigation  for 
many  years.  The  analytical  solutions  for  such  problems 
generally  include  three  basic  thermal  boundary  conditions: 

(a)  Uniform  heat  flux  per  unit  length  but  uniform  peri¬ 

pheral  wall  temperature  (Dirichlet  type  thermal 
boundary  condition) , 

(b)  Uniform  heat  flux,  axially  and  peripherally  (Neumann 
type  thermal  boundary  condition) , 

(c)  Uniform  wall  temperature  (Dirichlet  type  thermal 
boundary  condition) . 

It  has  been  shown  recently  that  the  point  matching 
technique  gives  an  approximate  solution  for  many  bound¬ 
ary-value  or  eigen- value  problems  involving  the  Laplace, 

Poisson  or  biharmonic  equations,  with  either  linear  or 

1* 


homogeneous  boundary  conditions.  Barta 


and  Conway 


applied  this  method  to  solve  a  number  of  torsion  and 
plate  problems  involving  triangular,  square,  hexagonal 


and  rhombic  shapes.  Sparrow 


has  accurately  determined 


the  flow  characteristics  for  laminar  flow  in  isosceles 


*  Numbers  in  brackets  denote  references  given  on  page  68  . 
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triangular  ducts  with  various  apex  angles.  Cheng 
applied  this  method  to  fully  developed  laminar  flow 
in  noncircular  passages  consisting  of  regular  polygons 
with  Dirichlet  type  boundary  conditions.  In  this  thesis, 
the  same  approach  is  utilized  for  fully  developed  lam¬ 
inar  flow  in  circular  ducts  with  diametrically  opposite 
flat  sides,  regular  polygonal  ducts  with  a  central  cir¬ 
cular  hole,  and  circular  ducts  with  a  central  regular 
polygonal  hole.  All  ducts  are  treated  for  Dirichlet- 
type  boundary  conditions  with  the  exception  of  regular 
polygonal  ducts  where  the  Neumann  boundary  conditions 
are  considered. 


The  variational  method  has  been  used  by  Sparrow  and 


Siegel 


to  solve  the  problem  of  rectangular  and  cir¬ 


cular-sector  ducts  with  various  thermal  boundary  conditions 
Because  of  its  success  in  the  above  problems,  the  method 
has  been  used  here  for  evaluating  the  heat  transfer  and 
flow  characteristics  for  circular  ducts  with  diametrically 
opposite  flat  sides  with  Dirichlet-type  boundary  conditions 
The  Raleigh-Ritz  variational  method  6  has  proved  to  be 
as  successful  as  the  point  matching  technique,  but  involves 
a  great  deal  of  work  if  accuracy  is  desired.  The  results 
obtained  using  this  method  are  found  to  be  quite  comparable 
with  those  computed  from  the  point  matching  technique. 

The  solutions  for  laminar  heat  transfer  for  both 
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parallel-plate  channels  and  circular  ducts  with  Dirichlet 
boundary  conditions  are  available  in  the  literature. 
Carter  7  and  Okubo  8  give  approximate  solutions  for 
the  Saint-Venant  torsion  problem  for  circular  cross 
sections  with  diametrically  opposite  flat  sides.  These 
approximate  results,  when  interpreted  in  terms  of  flow 
problems,  by  analogy,  substantially  agree  with  those  ob¬ 
tained  from  the  point  matching  and  variational  methods. 

A  comprehensive  study  on  heat  transfer  in  annular 
passages  for  hydrodynamically  fully-developed  laminar 
flow  with  arbitrarily  prescribed  wall  temperature  or 
heat  fluxes  is  presented  by  Lundberg,  McCuen  and 


Reynolds 


The  problem  of  bilateral  heat  transfer  in 


annuli  for  slug  and  laminar  flows  was  studied  by  Dwyer 


Gaydon  and  Nuttall 


11 


10 


established  the  upper  and  lower 
bounds  for  volume  flow  rate  through  cylindrical  tubes 
of  square  and  hexagonal  cross  sections  with  a  central 
circular  hole,  using  the  Schwartz  inequality  for  various 
hole  diameters.  The  present  study  for  doubly  connected 
boundaries  with  Dirichlet-type  boundary  conditions  is  based 
upon  the  method  of  point  matching.  Series  solutions 
have  been  utilized  by  satisfying  the  boundary  conditions 
at  the  inner  circular  boundary  exactly:  the  case  consid¬ 
ered  is  that  of  regular  polygonal  ducts  with  a  central 


regular  polygonal  hole 
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The  results  for  circular 
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ducts  with  a  central  regular  polygonal  hole  approximate 
very  closely  those  for  regular  polygonal  ducts  with  a 
central  circular  hole  when  the  number  of  sides  for  the 


4 


regular  polygon  approaches  a  very  large  number,  e.g. 


twenty. 


The  case  of  rectangular  ducts  with  uniform  peri- 
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in  the  form 


pheral  heat  flux  was  solved  by  Cheng 
of  an  infinite  series  using  separation  of  variables. 
Numerical  results  are  given  for  aspect  ratios  of  1,  2, 


4  and  10.  Sparrow  and  Siegel 


approached  the  same 


problem  with  aspect  ratios  of  1  and  10  using  the  varia¬ 


tional  method.  Eckert,  Irvine  and  Yen 
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presented 


the  heat  transfer  results  for  circular  sector  ducts 
with  apex  angles  ranging  from  11°  25'  20"  to  60°  using 
a  classical  analytical  method  for  both  Dirichlet  and 
Neumann  type  thermal  boundary  conditions.  Recently, 


Sparrow  and  Haji-Sheikh 
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presented  a  computation- 


oriented  method  of  analysis  for  ducts  of  arbitrary  shape 
with  arbitrary  thermal  boundary  conditions  by  using 
orthonormal  functions:  numerical  results  are  also  given 
for  ducts  of  circular-segment  cross  section.  The  method 
of  point  matching  has  also  been  used  for  evaluating  the 
heat  transfer  results  for  regular  polygonal  ducts  with 
Neumann  type  boundary  conditions  in  this  thesis. 
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CHAPTER  II 


THEORETICAL  ANALYSIS 


2 . 1  Governing  Equations 


The  governing  equations  and  the  associated  bound¬ 
ary  conditions  for  a  fully  developed,  constant  property 
and  non-dissipating  laminar  flow  are 
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(2) 


(a)  Dirichlet  type  boundary  conditions. 

u  =  9  =  0,  on  T 

(b)  Neumann  type  boundary  conditions. 


u  =  0 


dt 

dn 


~  ip,  on  r 


(3) 


(4) 


2  .  2  Series  Solution  for  Velocity 


The  general  solution  of  equation  (1)  in  polar  co¬ 
ordinates  (r,4>)  for  a  symmetrical  cross  section  (see 
Figures  1,  17,  32,  45)  is 
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u  = 


r  ’  2  A  dp' 
4  l|i  dz 


CO 


+  A  +  B  In  r 
o  o 


•♦I 


Ak  r  ' 


k  =1,2, . . . 


4*  r '  cos  k4> 


(5) 


The  above  solution  for  velocity  may  be  used  in  evaluating 
the  flow  characteristics  for  different  types  of  noncir¬ 
cular  ducts,  as  discussed  below,  with  no-slip  boundary 
conditions . 


(a)  Circular  Ducts  with  Diametrically  Opposite  Flat  Sides 


With  the  requirement  that  the  velocity  is  finite 
everywhere  and  in  particular  at  r '  =0  (see  Figure  1) 


it  follows  that 


From 


B  =  0, 
o  * 


du  _  n 
3$  "  0 


B,  =  0 
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at  *  -  f 


(6) 


k  =  2j  where  j  =  1,2,3,.,. 

Using  the  above  results  in  equation  (5) ,  the  velocity 
equation  is. 
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(b)  Regular  Polygonal  Ducts  with  a  Central  Circular  Hole 


The  velocity  distribution  for  regular  polygonal  ducts 
with  a  central  circular  hole  is 


u  = 


a2  dp 

d  dz. 
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o  o 


co 

I(s' 
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jN 


+  Bj  r  cos  jN<t> 


(8) 


The  unknown  coefficients  can  be  selected  so  as  to  satisfy 
the  boundary  condition  u  =  0,  on  the  inner  circular  bound¬ 
ary  (see  Figure  17) .  The  results  are 
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Substituting  Bq  and  B^  into  equation  (8) ,  we  get 
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(c)  Circular  Ducts  with  a  Central  Regular  Polygonal  Hole 


The  general  solution  of  equation  (1),  for  circular 
ducts  with  a  central  regular  polygonal  hole  is 
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To  satisfy  the  boundary  condition  u  =  0  on  the  outer 
circular  boundary  the  unknown  coefficients  must  be 


selected  as, 
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Substituting  the  above  constants  into  equation  (12) 
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(d)  Regular  Polygonal  Ducts 

If  the  velocity  is  finite  everywhere  and  particularly 
at  r  =  0  (see  Figure  45)  then  from  equation  (8) 
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and  the  remaining  equation  gives 
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2 . 3  Series  Solution  for  Temperature 


The  general  solution  to  the  temperature  equation  (2) 
can  be  expressed  as  the  sum  of  both  the  particular  and 
homogeneous  solutions 
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where  the  homogeneous  solution  of  the  temperature  equation 
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The  particular  solution  0^  can  be  obtained  by  selecting 
or  evaluating  some  function,  for  the  individual  case 
separately  as  discussed  below.. 


(a)  Circular  Ducts  with  Diametrically  Opposite  Flat  Sides 


With  the  velocity  field  given  by  equation  (7) ,  the 
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For  instance,  considering  only  the  jth  term  in  equation 
(7) ,  the  particular  solution  is 
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The  above  particular  solution  of  equation  (7)  and  the 
homogeneous  solution  represented  by  equation  (18) ,  pro¬ 
vide  us  with  the  temperature  distribution 
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(b)  Regular  Polygonal  Ducts  with  a  Central  Circular  Hole 


The  particular  solution  of  the  energy  equation  (2) 
can  be  found  with  the  velocity  field  given  by  equation 
(11) .  The  general  solution  for  equation  (2)  can  be  written 


as 
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On  substitution  of  and  B ^  into  equation  (22) ,  the 
temperature  difference  0  becomes, 
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’c)  Circular  Ducts  with  a  Central  Regular  Polygonal  Hole 


The  particular  solution  of  the  energy  equation  (2) 
can  be  found  by  using  equation  (14) .  Therefore  the  general 
solution  to  equation  (2)  can  be  written  as. 
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On  satisfying  the  boundary  condition  @  =  0  at  the  outer 
circular  boundary  and  setting  r’  =  c,  we  obtain, 
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On  substituting  A^  and  in  equation  (27) ,  the  temperature 
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difference  0  can  be  expressed  as 
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( d)  Regular  Polygonal  Ducts 


The  solution  to  equation  (2) ,  for  the  regular  polygonal 
ducts  can  be  expressed  as  follows 
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2 . 4  Analysis  by  the  Point  Matching  Method 


(a)  Evaluation  of  Unknown  Coefficients 


The  boundary  conditions  for  various  noncircular  ducts 
are  usually  difficult  to  satisfy  exactly.  However,  a 
completely  circular  boundary  provides  no  difficulty.  The 
other  boundary  e.g.  with  circular  ducts  with  diametrically 
opposite  flat  sides  and  regular  polygonal  ducts,  still 
requires  the  evaluation  of  the  unknown  coefficients  A ^ , 

B j ,  and  in  their  respective  equationsc 

In  solving  these  equations  one  might  apply  Fourier 


series 


Unfortunately  the  nature  of  these  equations 


precludes  the  use  of  Fourier  analysis  and  therefore 
some  other  approach  must  be  found  for  obtaining  these 


coefficients 


The  method  used  here  is  the  point  matching  method 
which  invokes  an  approximate  solution  by  taking  the  sum 
of  the  first  j  terms  of  an  infinite  series.  The  require¬ 
ment  that  the  boundary  condition  be  satisfied  at  discrete 
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points  on  the  remaining  unsatisfied  boundary  leads  to 

a  system  of  nonhomogeneous ,  algebraic  equations  for  the 

unknown  coefficients  A.,  B.,  A 1  and  B ,  which  can  be 

D  3  J  J 

easily  solved  using  a  digital  computer.  To  check  the 
accuracy  of  the  result  the  entire  procedure  may  be  re¬ 
peated  by  taking  more  points  at  the  boundary.  With  an 
increasing  number  of  points  at  the  boundary  one  can  ob¬ 
serve  the  results  approaching  a  definite  value  but  in 
practice  there  is  a  limitation  imposed  by  the  loss  of 
significant  figures  which  results  when  large  numbers  of 
linear  algebraic  equations  are  solved  numerically.  It 
is  also  worth  noting  that  the  values  of  these  coef f icients 
decrease  rapidly  with  an  increase  in  the  number  of  points 
taken  at  the  boundary. 

(b)  Unknown  Coefficients  for  the  Velocity  Equation 

An  approximate  solution  for  the  velocity  in  circular 
ducts  with  diametrically  opposite  flat  sides  is  obtained 
by  applying  equation  (7)  at  a  finite  number  of  points 
along  the  boundary  EFH  (see  Figure  1) .  The  appropriate 
number  of  points  on  the  straight  part  of  the  boundary  EP 
and  the  circular  part  FH  were  obtained  by  arbitrarily 
choosing  a  number  of  points,  on  the  circular  part  of  the 
boundary  and  varying  the  number  of  equiangular  spaced  points 
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on  the  straight  part  of  the  boundary,  or  vice-versa.  The 
numerical  results  for  velocity  and  shear  stress  for  ducts 


with  different  P's  were  evaluated  for  each  such  combina¬ 
tion  and  were  compared  with  the  results  of  Carter  and 
Oliphint  7  and  Okubo  8  .  The  boundary  errors  were 
also  calculated.  Table  1  represents  the  best  selection 
of  points  on  both  the  circular  and  straight  parts  of  the 
boundary .  Because  of  the  symmetry  of  the  cross  sections, 
equation  (7)  was  applied  to  only  one  quarter  of  the  cross- 
sectional  area  of  the  duct. 

For  regular  polygonal  ducts  with  a  central  circular 
hole,  the  boundary  condition  was  satisfied  exactly  at  the 
inner  circular  boundary.  The  boundary  condition  at  the 
remaining  regular  polygonal  boundary  was  satisfied  at 
some  equiangular  points.  For  this  shape,  because  of  the 
N-fold  symmetry,  only  one-2Nth  of  the  outer  boundary  need 
be  considered.  Experience  shows  that  a  relatively  greater 
boundary  error  exists  around  the  corner  point  and  that 
this  can  be  reduced  by  selecting  more  points  in  the 
vicinity  of  the  corner.  To  study  the  convergence  of  the 
solution  from  equation  (11),  4-point,  6-point,  8-point 
and  10-point  matchings  were  employed  for  various  values 
of  A  =  (a  cos  7r/N)/b. 


The  unknown  coefficients  for  equation  (14) ,  which 
represent  the  velocity  distribution  in  circular  ducts 
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with  a  central  regular  polygonal  hole,  can  be  evaluated 
by  the  point  matching  method  at  the  inner  regular  poly¬ 
gonal  boundary o  The  present  form  of  equation  (14)  is  based 
upon  the  condition  that  u  =  0  at  the  outer  circular  bound¬ 
ary  o  Due  to  the  symmetry  of  the  cross  section  of  the  duct 
only  one-2Nth  of  the  inner  boundary  was  considered*  To 
check  the  accuracy  of  the  results,  4-point  and  6-point 
matchings  were  employed  for  various  values  of  y  = 


*a  cos'f/N 

Since  the  computations  are  rather  extensive,  the  numerical 
results  for  the  unknown  coefficients  for  the  aforementioned 
cases  are  not  listed  in  this  thesis. 


Reference 


presents  the  numerical  values  of 


the  unknown  coefficients  A..  for  velocity  equation  (16) 
for  regular  polygonal  ducts.  These  values  were  also 
evaluated  with  equiangular  intervals  along  a  one-2Nth 


portion  of  the  boundary. 


Table  1  Number  of  Points  Matched  on 
Two  Boundaries  in  Circular  Ducts  with 

Diametrically  Opposite  Flat  Sides 
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[3° 

No.  of  Points  Matched  on 
Straiqht  Boundary 

No.  of  Points  Matched  on 
Curved  Boundary 

40 

17 

35 

50 

15 

31 

60 

17 

27 

70 

21 

23 

80 

30 

12 

85 

30 

6 

(b) 

Unknown  Coefficients  for  the 

Temperature  Equation 

Having  found  the  coefficients  and  for  their 
respective  equations  (7,11,14,16),  one  can  proceed  further 
to  evaluate  the  second  set  of  coefficients  A^  and 
from  the  temperature  equations  (21,26,29,30).  For  circular 
ducts  with  diametrically  opposite  flat  sides,  equation 
(21)  can  be  satisfied  at  the  same  points  as  selected 
in  evaluating  the  best  possible  values  for  the  unknown 
coefficients  of  the  velocity  equation  (7) .  For  regular 
polygonal  ducts  with  a  central  circular  hole  numerical 

results  are  obtained  for  the  unknown  coefficients  for 

% 

4-point,  6-point,  8-point  and  10-point  matching  for 
various  values  of  A  =  (a  cos  7r/N)/b,  by  satisfying  the 
temperature  equation  (26)  on  the  regular  polygonal 
boundary  of  the  duct.  The  boundary  condition  on  the 
circular  boundary  is  satisfied  exactly. 
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For  circular  ducts  with  a  concentric  regular  poly¬ 
gonal  hole  the  unknown  coefficients  and  A^  for  4-point 
and  6-point  matchings  are  obtained  for  various  values  of 
y  -  by  satisfying  equation  (29)  at  the  same 

points  on  the  inner  regular  polygonal  boundary  as  for 
velocity  equation  (14) .  Due  to  the  rather  extensive  num¬ 
erical  results,  the  values  of  the  unknown  coefficients  are 
not  presented  in  this  thesis. 

With  the  case  (b)  type  boundary  condition,  for  regu- 

/d©  '  \ 

lar  polygonal  ducts,  the  heat  flux  (— )  is  constant  on 
the  boundary.  This  makes  it  more  difficult  to  evaluate 
the  unknown  coefficients  for  temperature  equation  (30)  , 
The  following  procedure  was  utilized  to  evaluate  the  con¬ 
stant  Aj  for  different  sets  of  point  matchings  on  the 
regular  polygonal  boundary. 

The  temperature  gradient  normal  to  the  duct  boundary 
can  be  found  from 


<3© 

chT 


cos  <t>  -  (sin  4>/r)  for  (r/a)  -  (cos  tt/N)/cos  4>  0 


The  normal  temperature  gradient  is 
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A!  (jN)  h15"1  cos  (jN-l)t 


(31) 
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For  fully  developed  flow  with  uniform  heat  input  per 
unit  length,  the  axial  temperature  gradient  (dt/dz)  is 
constant.  Noting  the  following  relationships 


dt 
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it  can  be  shown  that 
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Substituting  equations  (32)  into  equation  (31) ,  the 
following  dimensionless  equation  can  be  written 
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,33) 


In  practice,  of  course,  the  finite  series  is  truncated 
and  the  boundary  condition  is  satisfied  at  a  finite  num¬ 
ber  of  suitably  selected  points  as  discussed  before c 
The  convergence  of  the  solution  can  be  demonstrated  by 
gradually  increasing  the  number  of  unknown  coef f icients „ 
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The  leading  term  of  the  homogeneous  solution  for  tempera¬ 
ture  5  A^  ( -  constant) ,  may  be  found  from  the  requirement 
that 

9'rdrd<J>  =  0  (34) 


Incidentally,  the  above  method  leads  readily  to  the  exact 
solutions  for  circular  ducts.  The  results  are, 
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The  unknown  coefficients  A^  for  various  regular  polygons 
were  obtained  by  numerical  integration.  For  a  square 
duct  the  result  checks  excellently  with  exact  integration 
using  Chebyshev  polynomials  for  cos  jN0  in  the  integration 
of  the  following  form, 


tt/N 


(cos  jN4>/cos^^  $) 


d<l> 


(37) 


The  solution  to  the  above  integral,  for  i^>  2,  is  pre¬ 
sented  in  the  appendix  whereas  for  i  =  2,  the  above  inte¬ 
gral  can  be  integrated  exactly.  The  formula  is 


cos  ^  <t>  cos  (j+l)^  d$ 


=  —  cos-1  0  sin  j4> 


(38) 


Tables  2  and  3  show  the  numerical  values  of  coefficients 
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Table 

2 

Coefficients 

A  . 
3 

for 

Regular  Polygonal 

Ducts 

2- 

a 

jN 

N 

Ao 

Ai 

2 

A3 

3 

0.833333  x 

10 

_1  -0.166667 

0.343873 

X 

io”15 

0.466933 

X 

io-15 
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0.147396 

-0.915117 

X 

10"1 

0.838082 

X 

10“2 

-0.271102 

X 

icf2 

5 

0.182372 
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X 

10"1 
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X 

io“2 

-0.312033 
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1<T2 
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io“2 

-0.226553 
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-2 

10  z 
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10 
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10 
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10”1 
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«2 

10  z 

-0.127912 
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10“2 

20 

0.245838 

-0.296296 
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10“2 

0.854921 
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ro 

1 
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tH 

-0.344408 

X 

io~3 

Table 

3  Coefficients 

A 

j  for 

Regular  Polygonal 

Ducts 
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A  * 
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Aj  and  A ^  for  regular  polygonal  ducts  of  N  =  3  to  10  in¬ 
clusive  and  N  =  20.  These  coefficients  were  obtained 
using  4-point  matching.  The  numerical  results  of  computa¬ 
tion  for  4-point,  8-point  and  12-point  matching  indicate 
that  the  leading  coefficients,  Aq  and  A^  for  velocity  and 
temperature,  respectively,  change  very  little.  This 
suggests  that  even  the  results  from  4-point  matching  are 
quite  close  to  the  exact  values.  These  leading  coefficients, 
for  equilateral  triangular  ducts,  do  not  change  at  all 
from  3  to  12-point  matching.  This  indicates  that  the 
3-point  matching  solution  for  equilateral  triangular  ducts 
is  very  close  to  the  exact  solution. 


2  o  5  Analysis  by  the  Variational  Method 


The  variational  method,  as  a  tool  for  the  solution 
of  laminar  heat  transfer  problems  with  various  boundary 
conditions,  is  discussed  fully  by  Sparrow  and  Siegel 
In  this  thesis,  only  the  circular  ducts  with  diametrically 
opposite  flat  sides  are  subjected  to  this  method  in 
order  to  evaluate  the  heat  transfer  characteristics 
using  case  (a)  type  boundary  conditions. 
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(a)  Velocity  Solution  for  Circular  Ducts  with  Diametrically 

Opposite  Flat  Sides 

In  solving  the  velocity  problem,  we  are  able  to 
benefit  from  the  knowledge  accumulated  in  the  theory  of 
elasticity.  There  is  a  direct  connection  between  equation 
(1)  and  the  differential  equation  for  the  stress  function 
in  the  Saint-Venant  torsion  problem.  Thus,  from  the  re¬ 


sults  given  by  Timoshenko 


one  can  write  a  first 


approximation  to  the  laminar  velocity  distribution  in 
the  circular  duct  with  diametrically  opposite  flat  sides 
(see  Figure  1)  as  follows, 


u 


2Pv 


'  2  ,  2\  /  2,2  2 

x  -  b  )  x  +  y  -  a 


(39) 


where  the  boundary  of  the  duct  is  given  by  the  equations 

x  =  +  b,  x^  +  y^  =  a^,  and  the  function  (x^  -  b^) 

2  2  2 

(x  +  y  -  a  )  assures  that  each  of  the  two  terms  satisfies 
the  condition  of  zero  velocity  at  the  boundary.  The  un¬ 
known  constant  Aq  is  to  be  determined  by  the  well-known 
Rit z  method  for  minimizing  a  variational  integral. 

Equation  (39)  is  introduced  into  the  following  variational 
integral 
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The  minimizing  process  is  carried  out  by  differentiating 
with  respect  to  the  unknown  constant  Aq  and  by  setting 
the  resulting  equation  to  zero,  obtaining  for  p  =  60°, 

Aq  =  (41) 

a 

To  test  the  accuracy  of  this  expression,  a  comparison  with 
the  velocity  solution  from  equation  (7)  is  made  in 
Figure  15  for  p  =  60° .  Profiles  are  plotted  along  the 
two  coordinate  axes.  The  agreement  is  generally  very 
good,  except  for  deviation  near  the  corner  and  the  centre. 

With  the  thought  of  improving  the  approximation, 
the  variational  method  was  re-applied  using  an  expression 
of  slightly  greater  complexity  than  equation  (39) „  The 
steps  in  the  calculation  will  be  shown  in  some  detail, 
so  that  future  calculations  which  proceed  in  the  same 
manner  may  be  described  only  in  outline  form.  The  pre- 
file  chosen  for  the  improved  approximation  is 

U  =  (”  2|i)  (x2”b2)  (x2+^2~a2)  (Ao  +  Alx2  +  A2y2) 

2  2  2  2  2 

The  factors  (x  ~b  )  and  (x  +y  -a  )  again  assure  that 
each  of  the  three  terms  satisfies  the  condition  of  zero 
velocity  at  the  boundary  described  by  x  =  +  b  and 
y  -  +  ja2-x2  0  The  constants  A  ,  ,  and  A^  remain  to 

be  determined. 

Substitution  of  equation  (42)  into  the  variational 
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integral  (40) ,  leads  to  a  result  which  can  be  further 
differentiated  with  respect,  to  each  of  the  A's  to  give 
us  the  values  of  Aq,  A^  and  A0  which  will  minimize  the 
value  of  I.  This  operation  yields. 


§£  =  °*  0.204  A  a2  +  0.00573  A, a4  +  0.03105  A  a4 

o  1  ^ 

+  0.206  -  0 

=  0;  0.00573  AQa2  +  0.008038  A1a4  -  0.004889  A2a4 

+  0.00939  =  0 

*  0;  0.03726  Apa2  +  0.000724  A^a4  +  0.02004  A2a4 

2 

+  0.03932  =  0  (43) 


Simultaneous  solution  of  this  group  of  linear  algebraic 
equations  provides  the  following  values  for  p  =  60° 


-  Q o 97466  =  OL 5527  =  0.13016 

2  s  Al  4  9  A2  4 

a  a  a 


(44) 


After  substituting  into  the  equation  (42) , 
u  =  (-  |i)(x2-b2)(x2+y2-a2^'0i97466  +  0.5527.x!  +  (45) 


The  mass  flow  rate  evaluated  from  equation  (45)  and  ex¬ 
pressed  in  terms  of  K  is  compared  with  that  obtained  from 
point  matching,  in  Figure  2.  The  agreement  is  generally 
very  good.  The  exact  solution  is  approached  more  and 
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mors  closely  as  more  terms  are  added  to  the  approximating 


polynomial,  for  u„  But  this  process  of  increasing  the 
number  of  terms  requires  a  tremendous  amount  of  time  to 
solve  the  variational  integral  (40)  as  the  terms  increase 


to  their  square  power „  The  other  improved  expression 
for  the  velocity  profile  is  as  follows. 


u 
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2HA  AX  y  ”a  A  o  1X  +A  y  +A3X  Y  +A4X  +A5y  ) 


(46) 


The  results  obtained  for  K  from  the  above  equation  are 
compared  in  Figure  2  with  the  value  of  K  obtained  from 
both  the  three  constants  and  point  matching .  The  agree¬ 
ment  is  found  to  be  very  good0 


(b)  Temperature  Solution  for  Circular  Ducts  with  Diametrically 

Opposite  Flat  Sides 


Consider  the  situation  of  uniform  temperature  around 
each  cross  section,  with  the  heat  transfer  per  unit  length 
in  the  axial  direction  constant .  Under  these  conditions, 
it  is  clear  that  the  wall  temperature  will  increase  from 
section  to  section  along  the  length  of  the  passage* 

For  this  situation  one  may  use  the  variational,  method 
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Table  4  Coefficients  A^  for  Circular  Ducts  with  Dia¬ 
metrically  Opposite  Flat  Sides-Variational  Method 
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80 

1.0211 

0.08094 

0.47216 

0.12757 

X 

102 

0.34725 

X 

10 

0.24920 

X 

10 

85 

1.0509  - 

•0.3078  - 

•0.33026 

0.16394 

X 

IQ2 

0.38020 

X 

10 

0.49344 

X 

10 

y  <, 

Table  5 

Coefficients  A j  for  Circular  Ducts  with  Diamet¬ 
rically  Opposite  Flat  Sides  -  Variational  Method 

(3° 

Ao" 

2 

a 

'■  if. . 

Ai  • 

4 

a 

A2  • 

4 

a 

A  »  . 

6 

a 

0 

< 

6 

a 

A5  ’ 

6 

a 

-1 

-1 

-1 

-1 

-1 

-3 

.5 

0.4660 

X 

10 

0.1462 

X 

10 

-0.1489 

X 

10 

-0.1477 

X 

10 

0.8328 

X 

10 

-0.5272 

X 

10 

-1 

-1 

-1 

-1 

-1 

-3 

10 

0.4745 

X 

10 

0.1487 

X 

10 

-0.1509 

X 

10 

-0.1393 

X 

10 

0.9406 

X 

10 

-0.6395 

X 

10 

-1 

-1 

-1 

-2 

-2 

20 

0.5045 

X 

10 

0.1647 

X 

10 

-0.1562 

X 

10 

-0.7097 

X 

10 

0.1123 

-0.1289 

X 

10 

-1 

-1 

-1 

-1 

-2 

30 

0.5370 

X 

10 

0.2026 

X 

10 

-0.1532 

X 

10 

0.1428 

X 

10 

0.1222 

-0.2900 

X 

10 

-1 

-1 

-1 

-1 

-2 

40 

0.5470 

X 

10 

0.2308 

X 

10 

-0.1257 

X 

10 

0.5023 

X 

10 

0.1088 

-0.5369 

X 

10 

-1 

-1 

-2 

-1 

-1 

-2 

50 

0.5083 

X 

10 

0.1857 

X 

10 

-0.6297 

X 

10 

0.8413 

X 

.10 

0.7611 

X 

10 

-0.7144 

X 

10 

-1 

-2 

-2 

-1 

-1 

-2 

60 

0.3990 

X 

10 

0.3475 

X 

10 

0.2945 

X 

10 

0.9350 

X 

10 

0.3228 

X 

10 

-0.8775 

X 

10 

-1 

-1 

-1 

-1 

-1 

-2 

70 

0.2267 

X 

10-0.1902 

X 

10 

0.1103 

X 

10 

0.6415 

X 

10-0.1571 

X 

10 

-0 . 6662 

X 

10 

-2 

-1 

-2 

-1 

-1 

-2 

30 

0.6264 

X 

10  -0.3619 

X 

10 

0.6371 

X 

10 

-0.1641 

X 

10 

-0.5899 

X 

10 

0.3619 

X 

10 

-2 

-1 

-3 

-1 

-1 

-2 

35 

0.1622 

X 

10 

-0.3762 

X 

10 

0.4506 

X 

10 

-0.7503 

X 

10 

-0.7490 

X 

10 

0.5397 

X 

10 
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Since  the  temperature  around  the  boundary  is  uniform, 
the  temperature  distribution  across  the  section  would 
be  expected  to  have  a  shape  not  too  different  from  that 
obtained  for  the  velocity.  Selecting  a  polynomial  for 
G  having  a  form  similar  to  equation  (46),  we  have, 


2 


9 


(48) 


Introducing  the  above  expression  for  9  into  the  varia¬ 
tional  integral  (47)  along  with  the  velocity  expression 

from  equation  (46)  and  carrying  out  the  integration  and 
minimizing  the  resultant  expression  by  differentiating 


linear  algebraic  equations  which  can  be  solved  to  obtain 

the  unknowns  A^,  A.J ,  AJ,  ...  A£  for  the  temperature  equation. 

The  unknown  coefficients  A  ,  A^,  A^  ...  A^  for  the  velocity 

equation  are  presented  in  Table  4  and  the  corresponding 

coefficients  A'  A'  A '  ...  A'  for  the  temperature  equation 

o J  1  ’  2  o 

are  presented  in  Table  5. 
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CHAPTER  III 


RESULTS  OF  ANALYSIS 


3 . 1  Circular  Ducts  with  Diametrically  Opposite  Flat  Sides 


(a)  Flow  Characteristics 


The  relationship  between  the  rate  of  flow  and  the 
pressure  drop,  which  in  dimensionless  terms  may  be  pre¬ 
sented  as  a  relationship  between  the  friction  factor  and 
Reynolds  number,  is  of  great  engineering  interest.  The 
mass  flow  rate,  w,  through  the  cross  section  of  the  duct 
can  be  calculated  from 

7 

o  o 


2  2 
a  -x 


u  dx  dy 


(49) 


This  can  be  further  subdivided  into  two  parts.  By  con¬ 
sidering  polar  coordinates  (r,<P)  the  mass  flow  rate,  w, 
becomes , 


w 


cos  {3 

cos  4> 


urdrd$ 


+ 


urdrd<t> 


7T 

2 


P 


o 


(50) 


or  w  =  pK  a^ ,  where  a  is  the  radius  of  the  cir¬ 

cumscribing  circle  and  K  is  a  constant.  The  values  of  K 
for  different  opening  angles  p  obtained  from  both  the  point 
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matching  and  variational  methods  are  tabulated  in 


Table  8,  together  with  results  for  rectangular  ducts 


(exact)  and  work  by  Carter 


.  Figure  2  represents 


a  comparison  of  the  constant  K  evaluated  by  the  varia¬ 
tional  method,  using  three  and  six  coefficients,  with 
values  obtained  by  point  matching. 

The  friction  f actor-Reynolds  number  relationship 
is  defined  by  the  following  equation 

i/i  ^ 

f  •  Re 


=  I  /  i 

2  v  10.  dzy 


u 


(51) 


m 


where  De  is  the  equivalent  hydraulic  diameter  and 
um  (*  w/pA)  is  the  average  velocity.  The  results  obtained 
from  equation  (51)  are  given  in  Table  9  and  are  graphi¬ 
cally  represented  in  Figure  14. 

For  Newtonian  fluids,  the  wall  shear  stress  is  de¬ 
fined  by  T  =  |i  (^f).  T'*ie  s^ear  stress  results  (from 
equation  (7))  for  the  straight  boundary  are  given  by 
the  equation 


T 


(-ac1) 


r  cos  $ 

2 


oo 

Y  A.  ( 2  j )  r2-^  1  cos  (2j-l)4> 


(52) 


and  for  the  curved  boundary  by 


oo 


j  =0,1,2 


2j-l 


cos  2 j 4> 


T  =  (-ac1) 


(53) 
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Table  6  Maximum  Velocity  for  Circular  Ducts  with  Dia¬ 


metrically  Opposite  Flat  Sides 


x  10 


p° 

Variational 

Method 

Point 

Matching 

Exact 

0 

2.5 

(circle) 

5 

2.508 

2.499 

— 

10 

2.502 

2.494 

— 

20 

2.463 

2.448 

— 

30 

2.340 

2.318 

— 

40 

2.082 

2.063 

— 

50 

1.673 

1.659 

— 

60 

1.139 

1.130 

1.139 

(rectangle) 

70 

0.577 

0.572 

0.573 

(  "  ) 

80 

0.1540 

0.1507 

0.1508 

(  "  ) 

85 

0.0399 

0.0380 

0.0380 

(  "  ) 

Table  7  Maximum  Shear  Stress  for  Circular  Ducts  with 
Diametrically  Opposite  Flat  Sides 


t/(  -  ac 


\ 


i/ 


x  10 


p° 

Variational 

Method  I- 

Point 

flat chine 

1  

Results  from 
Literature 

* 

Exact 

0 

5 

4.988 

CUR** 

5.259 

VER*** 

o 

•  i 

m 

(circle) 

10 

4.974 

II 

5.485 

II 

— 

— 

20 

4.889 

II 

5.808 

II 

— 

— 

28.95 

4.706 

II 

5.937 

If 

5.863 

- 

30 

4.677 

II 

5.941 

- 

— 

40 

5.742 

VER*** 

5.837 

II 

— 

— 

41.40 

5.788 

II 

5.802 

II 

5.888 

- 

50 

5.582 

II 

5.428 

II 

— 

- 

51.30 

5.492 

II 

5.348 

ii 

5.346 

— 

60 

4.631 

II 

4.624 

II 

4.650 

4.448 

( rectangle) 

70 

3.285 

II 

3.362 

It 

- 

3.226 

( rectangle) 

80 

1.732 

II 

1.732 

II 

- 

1 . 666 

(rectangle) 

85 

0.9071 

II 

0.8667 

II 

— 

0.7895 

(rectangle) 

*  Carter's  Results.  **  On  curved  boundary.  ***  On  straight  boundary. 
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Table  8  Constant  K  for  Circular  Ducts  with  Diametrically 
Opposite  Flat  Sides 


p° 

Variational 

Method 

Point 

Matching 

Results  from 
Literature 

*  Exact 

0 

tt/0.8 

(circle) 

5 

3.838 

3.925 

— 

— 

10 

3.818 

3.908 

- 

— 

20 

3.699 

3.773 

— 

— 

28.95 

3.424 

3.466 

3.393 

— 

30 

3.377 

3.416 

- 

- 

40 

2.779 

2.797 

- 

- 

41.40 

2.675 

2.691 

2.69 

- 

50 

1.963 

1.969 

— 

— 

51.30 

1.848 

1.854 

1.833 

— 

60 

1.098 

1.101 

1.095 

1.143 

(rectangle) 

70 

0.4119 

0.4125 

— 

0.4185 

(  ”  ) 

80 

0.06189 

0.06194 

— 

0.06217 

(  "  ) 

85 

0.00829 

0.00834 

0.00834 

(  "  ) 

Table  9  Values  of  f  •  Re  for  Circular  Ducts  with  Dia¬ 
metrically  Opposite  Flat  Sides 


p° 

Variational 

Method 

Point 

Matchinq 

Exact 

0 

mmt 

16 

(circle) 

5 

16.360 

15.998 

- 

10 

16.365 

15.980 

- 

20 

16.247 

15.933 

— 

30 

16.047 

15.862 

— 

40 

15.944 

15.842 

— 

50 

16.043 

15.986 

— 

60 

16.508 

16.475 

15.548 

(rectangle) 

70 

17.623 

17.603 

16.987 

(  "  ) 

80 

19.839 

19.822 

19.565 

(  "  ) 

85 

21.699 

21.551 

21.486 

(  "  .  ) 

*  Carter's  Results 
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The  results  from  the  variational  method,  for  the  curved 
and  the  straight  boundaries,  respectively,  using  equation 
(46)  are, 

1 

T  =  ci  (*2+y2)  ^x2-b2^Ao+A1x2+A2y2+A3x2y2+A4x4+A5y4^)  (54) 

and 

T  =  (bciXb2+y2-a2XAo+Aib2+A2y2+A3b2y2+A4b4+A5y4)  (55) 

The  shear  stress  distribution  along  both  the  straight  and 
the  curved  parts  of  the  boundary  is  presented  graphically 
in  Figures  3  and  4.  Table  7  gives  the  values  of  local 
maximum  shear  stress  on  the  boundary  for  different  opening 
angles  f3.  Figures  5,  6  and  7  represent  the  velocity 
distributions  along  three  different  directions  using  the 
point  matching  and  variational  methods  as  indicated  on 
the  respective  figures.  The  velocity  distribution  cal¬ 
culated  by  the  variational  method  is  smaller  in  magnitude 
than  the  point  matching  values  in  the  vicinity  of  the 
boundary  and  is  therefore  not  presented  in  Figures  6  and  7. 

(b)  Heat  Transfer  Characteristics 


With  the  velocity  and  temperature  fields  known, 
various  heat  transfer  characteristics  can  be  computed 
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readily.  The  limiting  Nusselt  number  (Nu)  based  on  the 


bulk  temperature  can  be  evaluated 


by. 


hD 

(Nu)lim  =  IT  =  -(D  /sa)  (bt/bz) 


ur ' dr ' d4> 


/ 


u9  r  '  dr  '  d<t> 


(56) 


The  temperature  distribution  using  the  point  matching 
method  in  all  three  directions  4>  =  0,  4>  =  0  and  4>  =  tt/2 
are  plotted  in  Figures  8,  9  and  10  respectively.  The  re¬ 
sults  using  the  variational  method  are  also  plotted  for 
comparison.  Figure  16  shows  a  comparison  of  the  tempera¬ 
ture  distributions  along  the  same  aforementioned  three 
directions  for  a  duct  with  an  opening  angle  of  (3  =  60° . 

The  use  of  the  variational  method  again  results  in 
only  minor  differences  (in  the  vicinity  of  the  duct 
boundary)  with  the  temperature  obtained  from  point  match¬ 
ing.  The  dimensionless  temperature  gradients  normal  to 
the  straight  part  of  the  boundary  obtained  from  equation  (21) 
can  be  expressed  as  follows, 
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and  from  equation  (48)  it  can  be  expressed  as, 


(■ 


b2-y2+a2  V  A^+A^b2+A2y2+A^b2y  2 


X- 


+A^b4+A£y4 


(58) 


The  above  results  for  various  values  of  (3  are  presented 
in  Figures  11  and  12.  The  results  from  the  variational 
method  generally  agree  with  those  of  the  point  matching 
method  except  for  ducts  with  (3  less  than  30  degrees.  This 
is  because  the  series  is  truncated  to  six  coefficients  in 
the  case  of  the  variational  method.  It  is  of  some  interest 
to  note  how  the  local  temperature  gradient  at  the  curved 
boundary  approaches  that  expected  for  a  circular  duct, 
as  (3  approaches  5  degrees. 

With  the  velocity  field  given  in  the  form  of  equations 
(7)  or  (46)  the  limiting  Nusselt  number  can  be  evaluated 
by  using  equation  (56) .  Exact  integration  can  be  performed 
using  equations  (46)  and  (48)  for  velocity  and  temperature 
respectively,  as  found  by  the  variational  method;  whereas 
in  equations  (7)  and  (21)  for  velocity  and  temperature 
as  evaluated  by  the  point  matching  method  it  was  found 


u0r '  dr'd<t>  is  very  difficult  to 


that  the  integral 


A 


evaluate  exactly  over  the  whole  area.  It  can  only  be 
applied  exactly  over  the  area  HOF  (see  Figure  1)  so 
numerical  integration  has  to  be  applied  for  the  rest  of 
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the  area,  EFO  (see  Figure  1) .  For  this  purpose  the  area 
EFO  of  the  cross  section  was  divided  into  square  grids 
with  width  (or  base)  and  height  no  greater  than  ^(1/50)  a 
cos  for  ducts  with  p  <^60°,  and  with  width  (or  base) 
and  height  no  greater  than  ^(1/100)  a  cos  p^  for  ducts 
with  (3  60°  . 


It  was  found  difficult  to  evaluate  the  limiting 
Nusselt  number  for  ducts  with  (3  =  85°  ,  using  the  point 
matching  method,  due  to  the  fact  that  the  boundary  error 
involved  was  much  higher  than  for  the  remaining  cases 
where  (3  <^85°.  The  resulting  Nusselt  numbers  using  the 
two  methods  are  shown  in  Figure  13  and  the  results  are 
listed  in  Table  11.  The  maximum  temperature  difference 
9  in  the  centre  of  the  cross  section  is  also  tabulated 
for  various  P's  in  Table  10. 

A  comparison  can  be  made  between  the  results  evaluated 
using  the  variational  method  and  the  results  evaluated 

using  the  point  matching  method.  In  the  case  of  the 

2  2  2  2  2 
variational  method,  the  functions  (x  -  b  )  and  (x  +  y  -  a  ) 

actually  satisfied  the  boundary  condition  u  =  0  =  0. 

Whereas  in  point  matching  there  always  exists  some  bound¬ 
ary  error  which  ultimately  effects  the  flow  and  heat 
transfer  results.  The  comparison  of  the  results  from  one, 
three  and  six  constant  A's  indicates  that  an  increase  in 
the  number  of  constants  using  the  variational  method  will 
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Table  10  Maximum  Temperature  for  Circular  Ducts  with 
Diametrically  Oppoiste  Flat  Sides 
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Table  11  Limiting  Nusselt  numbers  for  Circular  Ducts 
with  Diametrically  Opposite  Flat  Sides 
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not  improve  the  results  significantly,  compared  with 
the  amount  of  work  involved.  Hence  it  appears  that  the 
results  obtained  with  the  six  constants,  using  the 
variational  method,  are  better  than  the  results  obtained 
by  point  matching. 

3 . 2  Regular  Polygonal  Ducts  with  a  Central  Circular  Hole 

(a)  Flow  Characteristics 


The  mass  flow  rate  through  regular  polygonal  ducts 
with  a  central  circular  hole  can  be  calculated  from  the 
following  equation, 
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The  result  of  integration  after  substituting  the  velocity 
equation  (11)  is 
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In  (£2|_e) 

2A2  A 


I*  f  *)  { 


cos  (3  sin(jN-H)ft 

( jMH-1)  ( jlH-2) 


cos^^  (3  sin  i7r  ± 
(  j  N)  (  jN+  2)  Ajl^2 


n^N  sin(jN-l)(3  cos  (3 

( j  N- 1 )  ( j  N- 2 ) 


2jN 

n  sin  it r 


(jN)(jN~2)  cos^N  ^  p  A  jN+2 


(60) 


A  constant  K  is  defined  in  the  following  manner, 

/  4 


w 


=  pK  - 


M- 


dH 

dz 


The  value  of  K  for  various  values  of  A  were  determined 
by  matching  4,  6,  8  and  10  points  at  the  boundary.  On 
comparison  with  the  results  obtained  by  Gaydon  and  Nuttall 
11  for  the  upper  and  lower  bounds  of  K  for  square  and 
hexagonal  ducts  with  a  central  circular  hole,  it  was 
found  that  the  above  calculated  values  lie  in  between 
the  two  bounds.  The  results  of  10-point  matching  very 
closely  approximate  upper  bound  values  thus  demonstrating 
the  accuracy  of  the  results.  The  results  for  square  and 
hexagonal  ducts  together  with  those  for  upper  and  lower 
bounds  are  tabulated  in  Tables  12  and  13.  Also,  tabu¬ 
lated  in  Table  14  for  A  =  1.1,  20  and  1  to  10  inclusive 
are  the  results  for  the  constant  K  from  10-point  matching, 
for  various  regular  polygonal  ducts  with  a  central  circular 


(61) 
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Table  12  Values  of  the  Constant  K  for  a  Square  Duct 
with  a  Central  Circular  Hole 


A 

Lower 

Bound 

4-point 

6-point 

8-point 

10-point 

Upper 

Bound 

2 

0.02355 

0.02443 

0.02432 

0.02430 

0.02430 

0.02435 

3 

0.04322 

0.04423 

0.04410 

0.04408 

0.04408 

0.04410 

6 

0.06992 

0.07102 

0.07089 

0.07087 

0.07086 

0.07088 

10 

0.08320 

0.08434 

0.08421 

0.08419 

0.08419 

0.08420 

20 

0.09538 

0.09652 

0.09638 

0.09636 

0.09635 

0.09638 

GO 

0.13957 

— 

— 

— 

— 

0.14060 

- - - 

Table  13  Values  of  the  Constant  K  for  a  Hexagonal  Duct 
with  a  Central  Circular  Hole 


A 

Lower 

Bound 

4-point 

6-point 

8-point 

10-point 

Upper 

Bound 

2 

0.03566 

0.03709 

0.03693 

0.03690 

0.03689 

0.0369 

3 

0.07211 

0.07399 

0.07379 

0.07375 

0.07374 

0.07374 

6 

0.12313 

0.12544 

0.12518 

0.12513 

0.12511 

0  o 12516 

10 

0.14878 

0.15124 

0.15094 

0.15089 

0.15087 

0.15092 

20 

0.17224 

0.17480 

: 0.17448 

0.17442 

0.1744 

0.17443 

103 

0.21948 

0.2223 

0.2219 

0.2218 

0.2218 

0.22185 

106 

0.23788 

0.2408 

0.2404 

0.2403 

0.2403 

0.24036 

oo 

0.25622 

- 

— 

— 

— 

0.25892 

j  >il  »:tsu;>8  *.  yoi  M  imsianoD  drfi  ic  e^ricV  SI  eXcteT 

*-■  r  f.r.  )  J  ■  ,  rfiiw 


inxoq -01 

rnxoc-8 

Jmoc 

inxoq-- 

19WOJ 

i  nu( 

Ot  >£0 . 0 

0!  £S0.0 

S££so.o 

££££0.0 

eetso.o 

0X££0.0 

80  -£0 . 0 

0^  -0.0 

omo.o 

££££0.0 

£ 

86  VO.O 

0  ovo.o 

V80  0.0 

080 VO . 0 

S0XV0.0 

seeao.o 

a 

0£  £80 . 0 

0X£8O. 0 

x ;£tu.o 

£££80.0 

0££80.0 

ox 

88800.0 

e  -.660.0 

acaeo.o 

38)00.0 

sraec .0 

8£e80.0 

OS 

oacxx.o 

•- 

veexx.o 

ao 

h  s  'ic  -i  h  * r.. .  oD  o  -  9':  ■  -V  f  !.  elr:  - 

Xoli  f  roi  0  C  k  i  .:  leO  ri  i  r  w 


yeqqU 

bnpoB 

in  roq-8 

3n j  oq~d 

inioq-£ 

"ISWOvI 

bnuoS 

A 

98df 0 . 0 

oeaxo.o 

xeaxo.o 

eovf.o.o 

aaaeo.o 

S 

. 1 .0 

c :evo.o 

evcvo.o 

eeevo.o 

xxsvo.o 

axasx. o 

xxesx.o 

: IcSX.0 

exesx.o 

a 

eeoex.o 

£8i  ?  :.o 

£i.  et.o 

6Y8£X . 0 

ox 

£££V  X . 0 

£M  1.0 

8££V I .0 

08£V  X . 0 

£££YX . 0 

os 

8XSS.0 

exss.o 

£SSS .  0 

3£0X£ . 0 

£0££ . 0 

£0  .0 

80 i  S  .  0 

88Y££ . 0 

°ox 

- 

— 

Sf  o r  S  .0 

oo 

42 


hole.  The  graphical  results  are  shown  in  Figure  18. 

The  results  for  K  using  A  =  1  cannot  be  obtained 
from  equation  (11)  since  the  solution  is  only  good  for  a 
geometry  with  N-fold  symmetry.  However  the  case  of  A  =  1.1 
was  considered. 


With  the  velocity  distributions  represented  by 
equation  (11)  in  the  form  of  finite  series,  wall  shear 
stress  distributions  can  be  computed  readily.  In  the 
case  of  Newtonian  fluids  the  shear  stress  is  calculated 
from  t  =  ix(du/dn)  and  du/dn  =  (sin  <t>)/r  •  du/d4>  -  cos  $ 
du/dr,  which  in  turn  gives 


T 


2  \ 
— a - ) 

4r  In  t]  / 


cos 


<t>  + 


Aq  cos  $ 
r  In 


CO 

-  V  Aj  (jN) 
3 “1,2,3 


r-^N  ^  cos  (jN~l)0 


oo 


l 


A^jN)  T!2jN 
~  jN+1 


cos  ( jN  +  1)$ 


(62) 


Shear  stress  distributions  at  the  outer  boundary  for 
N  =  4  are  shown  graphically  in  Figure  19.  It  was  found 
that  the  shear  stress  distributions  along  the  inner 
circular  boundary  are  fairly  uniform  for  most  of  the 
cases  and  are  therefore  not  presented. 

A  friction  f actor-Reynolds  number  relationship  can 
also  be  defined  for  regular  polygonal  ducts  with  a  central 
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Table  15  Values  of  f  •  Re,  for  Regular  Polygonal  Ducts  with 
Central  Circular  Hole,  10-point  Matching 
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circular  hole  by  equation  (51) .  The  friction  factor- 
Reynolds  number,  f  •  Re,  is  consistent  with  that  for  a 
simply  connected  duct;  however  it  cannot  be  related  to 
the  wall  shear  stresses  in  a  multiply  connected  duct. 

It  is  considered  to  be  more  useful  if  the  friction  factor 
for  each  closed  independent  boundary  is  defined  by  finding 
the  location  of  maximum  velocity  from  equation  (11)  for 
a  given  angular  position c  But  this  location  cannot  be 
found  precisely  and  hence  the  job  of  obtaining  a  friction 
factor  for  each  boundary  is  a  difficult  one  and  is  not 
included  in  this  study e  The  numerical  values  of  f  •  Re 
from  10-point  matching  are  listed  in  Table  15  for  various 
cases.  These  graphical  results  are  also  shown  in 
Figure  20 » 

(b)  Heat  Transfer  Characteristics 

The  temperature  gradient  at  the  wall  is  of  general 
practical  interest.  By  considering  the  case  (a)  type 
boundary  condition ,  the  temperature  gradient  along  the 
inner  and  outer  boundary  can  be  evaluated  easily.  The 
temperature  gradient  normal  to  the  regular  polygonal 
boundary  can  be  expressed  .by  using  equation  .26;  as 
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The  temperature  gradient  for  various  values  of  A  along 
the  regular  polygonal  boundary  are  presented  in  Figure  21 
for  the  case  N  “  4  only.  The  computed  value  for  the 
temperature  gradients  at  the  circular  boundary  are  found 
to  be  fairly  uniform  and  are  therefore  not  presented. 

For  triangle  and  square  the  temperature  gradient  at  the 
corner  is  zero.  It  is  noted  that  the  convergence  of  the 
heat  transfer  characteristics  to  the  simply  connected 
case  (A  *  oo)  is  very  slow. 

Nusselt  numbers  for  regular  polygonal  ducts  with 
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a  central  circular  hole  can  be  evaluated  by  using  the 
velocity  and  temperature  equations  (11)  and  (26)  . 

The  temperature  and  velocity  distributions  are 
presented  graphically  in  Figures  22  to  29  inclusive  by 
using  10-point  matching  along  the  two  directions  0  =  tt/N,  and  0 
for  square  and  20-sided  ducts. 

Comparison  is  also  made  with  the  results  from 
regular  polygonal  ducts  with  no  central  hole.  The  re¬ 
sults  from  concentric  annuli  are  presented  in  Figure  23, 

The  velocity  distribution  along  0  =  tt/20  for  N  =  20  and 
A  =  2  is  compared  with  the  corresponding  one  for  an  annulus 
in  Figure  30 „  The  agreement  is  quite  good  up  to  (r°  -  b) / 

(a  -  b)  =0.35;  but  from  there  on  the  velocity  for  the 
20-sided  regular  polygon  is  slightly  lower  than  for  the 
annulus.  This  observation  could  have  been  anticipated, 
since  one  would  not  expect  any  effect  from  an  outside 
regular  polygonal  boundary  in  the  vicinity  of  an  inner 
circular  boundary  for  a  reasonably  large  value  of  A. 

A  heat  transfer  coefficient  can  be  defined  at  each 
boundary  by  using  the  fluid  bulk  temperature  of  the 
flowing  fluid  between  the  respective  boundary  and  a 
closed  line  of  minimum  temperature.  By  the  trial  and 
error  method  the  location  of  the  minimum  temperature 
can  be  found  from  dQ/dr“  =  0  using  equation  (26)  for  a 
particular  angular  position.  Numerical  integration  will 
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be  required  to  evaluate  the  Nusselt  number  for  each  in¬ 
dependent  boundary o 

An  overall  average  coefficient  may  be  defined  in 
terms  of  the  uniform  wall  temperature  and  the  bulk  fluid 
temperature  for  the  total  heat  transfer  rate  from  both 
walls.  The  Nusselt  number  can  be  evaluated  in  a  manner 
similar  to  that  used  for  the  simply  connected  ducts  4 
using  the  thermal  boundary  conditions  presently  under 
consideration „  Therefore  equation  (56)  is  applied  in  this 
case  to  evaluate  the  limiting  Nusselt  number  for  various 


values  of  A.  The  integral  J  J  ur'dr‘d4>  was  solved 
analytically,  whereas  the  solution  for  the  integral 

u0r'dr8d4>  was  found  by  numerical  means.  The  results 


obtained  from  8-point  matching  together  with  results  for 


an  annulus 


10 


are  plotted  in  Figure  31  and  the  results 
obtained  by  the  point  matching  method  are  tabulated  in 
Table  16 „  Figure  31  reveals  that  the  Nusselt  number 


curve  approaches  that  obtained  by  Dwyer 


10 


for  higher 


values  of  A.  It  is  to  be  expected  that  the  Nusselt 
number  for  the  20-sided  regular  polygonal  duct  with 
a  central  circular  hole  and  A  =  oo  would  be  comparable 
to  the  Nusselt  number  of  4.36  obtained  for  a  circular 
duct  by  classical  analytical  methods.  Nusselt  numbers 
are  also  calculated  for  N  =  4  for  very  large  values  of 


A.  For  A  =  1000  and  A  =  10  ,  the  Nusselt  numbers  are 
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found  to  be  4.556  and  4.056  respectively.  The  Nusself 
number  for  N  =  4  with  no  inside  circular  hole  is  3o6089 
which  indicates  that  the  above  convergence  is  very  slow. 

3 . 3  Circular  Ducts  with  a  Central  Regular  Polygonal  Hole 

(a)  Flow  Characteristics 


The  mass  flow  rate  can  be  calculated  by  using  an 
equation  similar  to  equation  (59) ,  which  is 
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After  substitution  of  u  from  equation  (14)  and  on  inte¬ 
grating,  we  get 
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Defining  a  constant  L  in  the  same  manner  as  in  equation 
(61)  we  get, 

-  -  K-t  S) 

The  results  obtained  for  L  from  the  above  equation  for 


,  closely  approximate 


N  =  20  with  various  y  =  - 

those  obtained  for  concentric  cylinders „  This  serves  to 
demonstrate  the  accuracy  of  the  results „  The  numerical 
results  for  constant  L  with  various  values  of  y  using 
4-point  and  6-point  matchings  for  various  circular  ducts 
with  a  central  regular  polygonal  hole  are  given  in 
Tables  17  and  18^  the  graphical  results  for  6-point 
matching  are  plotted  in  Figure  33. 

The  shear  stress  at  the  inner  regular  polygonal 
boundary  can  be  represented  using  equation  (14) , 
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Graphical  results  for  the  shear  stress  distribution  at 
the  inner  boundary  are  shown  for  y  -  2  only,  in  Figure  34 a 
The  graphs  show  a  large  fluctuation  which  might  be  sus¬ 
pected  due  to  the  large  boundary  error  at  the  inner  regu¬ 
lar  polygonal  boundary  for  lower  values  of  N.  With  the 
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Table  17  Values  of  the  Constant  L  for  Circular  Ducts  with 
Central  Regular  Polygonal  Hole,  4-point  Matching 
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Table  18  Values  of  the  Constant  L  for  Circular  Ducts  with  a  Central 
Regular  Polygonal  Hole,  6-point  Matching 
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increase  in  the  number  of  points  at  the  boundary,  the 
unknown  coefficients  for  the  velocity  equation  start 
diverging  in  absolute  value  which  resulted  in  tremendously 
high  values  for  the  velocity  distribution.  The  reason 
is  unknown  to  the  author. 

It  was  found  from  equation  (67)  that  the  shear 
stress  distributions  along  the  outer  boundary  are  fairly 
uniform  for  most  of  the  cases  and  therefore  not  graphi¬ 
cally  presented. 

Equation  (51)  was  used  to  define  the  friction 
f actor-Reynolds  number  relationship.  In  this  case  too 
it  was  physically  more  meaningful  to  define  one  friction 
factor  for  each  independent  boundary.  But  with  the 
velocity  field  given  by  equation  (14)  it  was  found  very 
difficult  to  establish  a  position  for  the  maximum  vel¬ 
ocity  for  a  particular  angular  position.  Therefore 
no  study  was  done  for  this  type  of  friction  factor.  The 
numerical  results  for  the  friction  f actor-Reynolds  num¬ 
ber  were  computed  for  various  circular  ducts  with  a  cen¬ 
tral  regular  polygonal  hole  for  y  =  1  to  9  inclusive 
and  are  listed  in  Table  19,  for  6-point  matching.  The 
graphical  results  are  also  presented  in  Figure  35. 

(b)  Heat  Transfer  Characteristics 


The  temperature  gradient  normal  to  the  inner  boundary, 


edi  ^xsbnnod  9di  2b  ain xoq  2o  leaiinun  sdi  nx  eassioai 

^TEte  noxi  >t/p©  y^xooIqv  ©rf  xo2  ain©x:>I22©oo  nwondnu 
Ylauobfi  iwaiJ  nx  beiljaex  rfDi.n  *  suJLjsv  ait/Ioads  nx  pnipxsvib 
.ic  691  9  :T  .noxindxxiax >  vjxo of  >v  ait  ,•  10^  asblxv  rfpxd 

.xodii/s  ©di  oi  owcnrfnu  ax 
xx  ada  ©di  xdi  (Vt  iOxisups  n:oxi  bnoo2  xw  21 
Y1xxx2  ©xe  '  xebnucc.  isjuo  edi  pnols  anoidjjdxxiaxb  aaexiia 
-j  rfq  :ip  ion  ©xo2©xed  r  bns  a  axo  ©di  2o  iaora  xo:;  mxoixcu 

.beinaeoxq  \  Leo 

noxioxxi  adi  oniis^  o 2  beau  "xw  (Ic)  naxixnpa 
ooi  98£o  axdi  nl  .qxrfenoxixlsx  xedmun  ablonys/f  iotOBr; 
no.  ioi  3  eno  9n  2 3b  ‘bn  r '?9fn  e-  . im  Y  snxeYdq  f  w 

edi  dixw  ij?e  . Yxermifod  o  t  bnaqsbni  r.->£9  xo  xoxos2 
vxsv  bru/o2  esw  2x  ( ^  I )  noi  Gt/p©  Y<*  fl  >vxp  b  9x2  y^XdoXdv 
-lev  miJ/nxxsm  edi  xc  t  no  ..  x  :o;  i  ra  :L(  &  £>  oi  It/c  22xb 
9Xo29X©dT  .noxiiaoq  xslupn  ix.  no  t  ir  -q  x  xo:  Y  "Xoo 
.  xoiox2  noxioxxi  2v  eqYi  aidi  xo  nob  ssw  y^^8  orI 
-mc/n  ablon  9fl-~oiox2  n  on  ©r  *  i o:  jiluaax  f  >  oxxe/iux 
-neo  s  dixw  aioo  >  xxlnoxxo  ano.  xxv  -io'  £>b2u<  ioo  ©  aw  X9d 
9vis»iXonx  6  oi  X  =*  y  bXori  Ism  pvloq  xxXupex  Xsxi 
sriT  .pnxrfoixm  Inxoq  d  xcl  *£1  ©lix  .  nl  b  riei  6l«  bnx 
,d£  euoi  l  t  .c  by  .  lea::  ;  oa  it  •.  ■  c/a©x  Isoxrfqs.p 

aoiOa -,is  ox  xdP  x ».ir  n x xT  iseH 

.vxxbnjod.  xonni  9rii  oi  Xsmxon  in9xbxxp  ©XLisxeqmei  ©dT 


56 


with  the  thermal  boundary  condition  of  circumferentially 
uniform  wall  temperature  at  a  given  axial  position  and 
axially  uniform  but  unequal  heat  fluxes  from  the  inner 
and  outer  walls,  can  be  obtained  by  using  equation  (29) , 
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The  results  from  the  above  equation  for  the  case  of 
7=2  only 9  for  various  regular  polygonal  boundaries,  are 
presented  in  Figure  36.  The  large  fluctuation  in  the 
temperature  gradient  was  expected  due  to  a  limited  number 
of  points  matched  at  the  inner  regular  polygonal  boundary 
and  also  the  effects  from  the  large  errors  in  velocity. 


(68) 
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It  was  found  that  the  temperature  gradient  for  the  circular 
boundary  was  fairly  uniform  for  most  of  the  cases.  The 
velocity  and  temperature  profiles  along  the  two  directions 
0  =  0  and  0  —  tt/N  are  shown  in  Figures  37  to  40  inclusive 
for  N  =  6  and  9  for  various  values  of  yc  Figures  41  and 
42  represent  the  velocity  distribution  along  0  =  0  and  0  =  7r/10, 
respectively  (for  N  =  10  only)  for  various  values  of  ya 
A  comparison  is  also  made  with  the  results  obtained  for 
the  circular  cross-section ,  The  temperature  and  velocity 
profiles  in  the  other  directions  are  bounded  by  the  two 
profiles  for  0  =  0  and  0  =  tt/N,  The  velocity  distribution 
along  0  =  tt/20  for  N  =  20  and  7  =  2  is  compared  with  the 
corresponding  one  for  an  annulus  in  Figure  43 « 

With  the  case  (a)  type  boundary  condition  and  an 
average  heat  transfer  coefficient  ,  which  could  be  defined 
at  each  boundary  by  using  the  bulk  fluid  temperature  of 
the  fluid  flowing  between  the  respective  boundary  and  a 
closed  line  of  minimum  temperature ,  the  Nusselt  number 
can  be  expressed  by  equation  (56) 0  It  could  be  expected 
that  numerical  integration  would  be  required  for  the 
calculation  of  the  Nusselt  number  for  each  independent 
boundary o 

The  Nusselt  number  results  which  were  computed  by 
numerical  integration  for  y  =  2  to  9  inclusive  and  for 
N  =  4  to  9  inclusive  are  tabulated  in  Table  20  and  some 
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plotted  in  Figure  44 „  These  results  seem  to  be  in  good 


agreement  with  those  of  Dwyer 
values  of  7  and  Nc 


10 


for  annuli  with  large 


3  o  4  Regular  Polygonal  Ducts 


(a)  Flow  Characteristics 


The  numerical  results  for  fully  developed  .laminar 
flow  in  regular  polygonal  ducts  have  been  discussed 


by  Cheng 


using  a  series  solution,  so  these  will 


not  be  discussed  here. 


(b)  Heat  Transfer  Characteristics 


Using  the  temperature  distribution  defined  by 
equation  (16) ,  the  Nusselt  number,  average  wall  tempera¬ 
ture,  and  bulk  fluid  temperature  may  be  computed  readily 
and  are  defined  as  follows. 
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The  calculated  results  using  12-point  matching  for  the 
average  wall  temperature,  bulk  fluid  temperature,  and 
normal  temperature  gradient  at  the  wall  are  listed  in 
Table  21  for  various  polygons «  To  check  the  convergence 
of  the  solution,  4-point,  8-point  and  12-point  matchings 
were  employed  for  all  computations  involving  the  above 
heat  transfer  characteristics „  The  numerical  results 
show  that  the  differences  among  the  above  three  kinds 
of  matching  are  negligible;  for  example,  the  temperature 
gradients  at  the  wall  generally  agree  to  four  significant 
figures 0  This  indicates  that  the  boundary  errors  are 
very  small  and  that  the  solution  converges  very  quickly . 

The  results  for  the  Nusselt  number  and  average  wall  temp¬ 
erature  are  tabulated  in  Tables  22  and  23  respectively. 

The  numerical  values  for  fluid  bulk  temperature  show 
nearly  the  same  trend  as  the  values  for  average  wall 
temperatures,  as  shown  in  Table  23 .  Temperature  distri¬ 
butions  from  12-point  matching  in  two  directions  <P  ®  tt/N  and 
-0  are  plotted  in  Figures  47  and  48,  respectively, 
for  several  representative  polygons . 

The  question  arises  as  to  the  accuracy  of  the  Nusselt 
numbers  calculated  for  different  regular  polygonal  ducts. 

A  comparison  can  be  made  between  the  results  calculated 
by  point  matching  and  the  results  available  in  the  litera¬ 
ture  »  For  instance  the  Nusselt  number  obtained  for  a 
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Table  21  Heat  Transfer  Results  for  Regular  Polygonal 
Ducts  with  Uniform  Peripheral  Heat  Flux 
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20-sided  regular  polygonal  duct  can  be  compared  with 
that  of  a  circle c.  An  approximate  0o8  per  cent  difference 
between  the  two  values  indicates  that  they  are  in  close 
agreement o  Also,  the  Nusselt  number  obtained  for  a 
square  duct  is  in  close  agreement  with  values  calculated 


by  Cheng 


13 


and  Sparrow  and  Siegel 


The  square  duct  with  uniform  peripheral  heat  flux 


was  studied  by  Cheng 


13 


using  separation  of  variables 


and  this  led  to  rather  involved  results  in  the  form  of 
an  infinite  series.  The  Nusselt  number  of  3.82,  com¬ 
puted  by  him,  appears  inconsistent  with  his  results  for 
rectangular  ducts  with  a  higher  aspect  ratio,  and  hence 
may  be  in  error.  The  result  from  this  investigation, 
3.091,  checks  exactly  with  Sparrow  and  Siegel ’ s  5 
estimation,  using  the  variational  method,  to  three  sig¬ 
nificant  figures o  All  values  for  the  Nusselt  number  for 
different  regular  polygonal  ducts  lie  between  10892  and 
4.364,  where  the  first  value  corresponds  to  an  equilateral 
triangular  duct  and  the  latter  value  to  a  circular  duct. 
The  Nusselt  number  of  1.892  can  be  compared  with  that 


of  2 o 3691  as  computed  by  Eckert,  Irvine  and  Yen 


14 


for  a  circular  sector  duct  with  an  apex  angle  of  60  de= 


grees  under  case  (b)  type  boundary  conditions.  Here 
the  difference  between  Nusselt  numbers  is  evident  and 


is  due  to  the  circular  boundary  involved.  The  average 
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Table  23 


Average  Wall  Temperature  for  Regular  Poly¬ 
gonal  Ducts  with  Uniform  Peripheral  Heat  Flux 
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Nusselt  numbers  for  uniform  peripheral  wall  temperature 
and  uniform  peripheral  heat  flux  are  also  presented  in 
Figure  46  for  comparison. 

The  locations  of  the  "hot  spots"  in  equilateral 
triangular  ducts  are  indicated  in  Figure  49,  which  re¬ 
presents  the  temperature  distribution  along  the  sides  of 
different  regular  polygonal  ducts.  The  temperature  dif¬ 
ference,  ©^  is  negative  at  the  centers  of  the  sides.  This 
simply  means  that  the  temperatures  at  these  points  are 
lower  than  the  average  temperature,  t^,  of  the  fluid.  These 
effects  are  quite  important  in  thermal  stress  analysis. 
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3 . 5  Concluding  Remarks 
General  Remarks 

1.  The  results  of  this  analysis  will  apply  only  to  fully 
developed  flow.  That  is,  the  asymptotic  forms  of  the 
Navier-Stokes  and  energy  equations  given  by  equations 

(1)  and  (2)  are  assumed  valid.  This  will  apply  at  points 
well  downstream  of  the  entrance. 

2.  The  "hot  spots"  at  the  corners  of  the  regular  polygonal 
ducts  with  uniform  peripheral  heat  flux  are  primarily 
caused  by  the  rather  low  velocity  field  existing  near 
the  corner.  In  reality  heat  conduction  along  the  bound¬ 
ary  exists  and  this  tends  to  decrease  the  temperature 

at  the  corners. 

3.  The  results  of  computation  show  that  higher  boundary 
errors  exist  at  the  corners  of  the  ducts.  The  exact 
reason  for  this  is  not  known  to  the  author.  However 
these  effects  are  limited  to  the  corners  only  and  have 
negligible  effect  on  the  overall  flow  or  heat  transfer 
characteristics . 

Detailed  Remarks 

1.  The  present  study  indicates  that  both  the  point 

matching  and  variational  methods,  for  the  approximate 

solution  of  the  stated  boundary  value  problem,  appear  to 
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have  advantages  over  the  finite  difference  technique  for 
the  problems  considered.  In  both  the  point  matching  and 
variational  methods  the  accuracy  of  the  result  was  found 
to  depend  upon  the  number  of  unknown  coefficients  selected. 

Other  techniques  of  satisfying  boundary  conditions  such 
as  least  square  method  are  not  considered  in  this  thesis. 

2.  A  computer  programme  can  be  easily  prepared  to 
obtain  results  using  the  point  matching  technique,  whereas 
with  the  variational  method  a  great  deal  of  work  is  in¬ 
volved  in  attaining  the  same  degree  of  accuracy. 

3.  With  the  velocity  field  given  in  the  form  of 
series  solutions  (equation  5) ,  the  particular  solution 
for  the  temperature  can  always  be  found  using  the  point 
matching  technique.  The  temperature  distribution,  using 
the  variational  method  (for  circular  ducts  with  diametri¬ 
cally  opposite  flat  sides) ,  can  also  be  defined  in  the 
same  form  as  that  of  the  velocity  distribution. 

4.  The  point  matching  method  can  be  used  for  any 
compact  section  but  elongated  sections  may  give  rise  to 
divergent  solutions  or  physically  unreasonable  results. 

5.  For  geometrical  configurations  with  both  curvilinear 
and  straight  line  parts  forming  the  boundary,  it  was  very 
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difficult  to  select  the  appropriate  number  of  points 
on  both  boundaries  to  give  the  best  possible  heat  transfer 
results  coincident  with  a  fairly  small  boundary  error, 

6.  The  point  matching  method  can  also  be  applied  to 
regular  polygonal  ducts  for  fully  developed  laminar  heat 
transfer  with  Neumann  boundary  conditions,,  The  numerical 
results  were  found  to  converge  quickly  to  the  exact  values „ 
It  is  believed  that  the  Nusselt  number  results  obtained 
for  this  case  are  well  within  0,5  percent  of  the  exact 
values . 

7,  The  question  arises  as  to  the  accuracy  of  the  vari¬ 
ous  calculated  results.  The  results  for  the  mass  flow 
rate  shown  in  various  tables  are  very  close  to  the  exact 
values.  For  example,  the  mass  flow  rate  values  for  a 
20” sided  regular  polygonal  duct  with  a  central  circular 
hole  or  a  circular  duct  with  a  central  regular  polygonal 
hole  approach  that  of  an  annulus  for  the  same  values  of 
A  or  y„  For  circular  ducts  with  a  central  regular  poly¬ 
gonal  hole,  the  calculated  values  for  the  friction  factor 
and  the  shear  stress  seem  to  differ  significantly  from 
the  actual  values,  the  reason  being  large  errors  in 
velocity  and  temperature  at  the  boundary.  The  boundary 
errors  calculated  for  velocity  and  temperature  for  cir¬ 
cular  ducts  with  diametrically  opposite  flat  sides,  regu¬ 
lar  polygonal  ducts  with  a  central  circular  hole  and 
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regular  polygonal  ducts  are  found  to  be  very  small „  The 
maximum  error  is  estimated  to  be  less  than  one  percent 
of  the  maximum  value  of  velocity  or  temperature  in  the 
duct  „ 


8.  It  appears  possible  to  extend  the  present  method 
of  point  matching  to  cases  with  uniform  heat  sources  and 


viscous-dissipation  effects 


17 


whereas  the  variational  method 


cumbersome . 


The  variational  method 


for  the  shapes  discussed. 


would  be  extremely 


can  be  applied  to 


laminar  flow  problems  in  circular  ducts  with  diametrically 
opposite  flat  sides  with  case  (a)  and  (b)  type  boundary 
conditions . 

10.  The  point  matching  technique  may  also  be  applied 
to  the  four  geometrical  shapes  discussed  with  various 
other  combinations  of  thermal  boundary  conditions 
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APPENDIX 


Exact  Integration  using  Chebyshev  Polynomials 

As  there  is  no  exact  integration  available  for 
the  following  integral 


o 


where  i  ^>  2 ,  it  was  difficult  to  evaluate  the  unknown 


coefficient  Aq  by  equation  (34)  for  regular  polygonal 


ducts  with  the  case  (b)  type  boundary  condition.  The 
solution  to  the  above  integral  was  attained  using 
Chebyshev  polynomials  which  changed  the  above  integral 
into  a  simpler  form,  the  solution  to  which  exists  in 
the  literature.  This  simplified  integral  could  then  b 
evaluated  using  a  digital  computer. 

The  following  conversion  was  adopted  to  evaluate 
the  simpler  form  of  the  integral, 


(1) 


Nj/2 


(-1) m(Nj-m-l) l 

ml  (Nj-2m) 1 


m  =0 


(2) 


Substituting  in  the  above  integral,  one  gets 


2N j-2m 


(3) 


m  =  0 


o 


which  is  the  required  form  to  be  integrated  exactly,  the 
result  evaluated  via  the  digital  computer. 
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Fig.  1  Coordinate  system  for  circular  duct 

with  diametrically  opposite  flat  sides. 
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Fig.  2  Mass  flow  rate  for  circular  ducts  with  diametrically  opposite 
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Fig*  10  Temperature  profiles  along  <t>  =  tt/2  in  circular  ducts  with  diametrically 
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Fig.  17  Coordinate  system  for  N-sided  regular  polygonal 


duct  with  a  central  circular  hole. 
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Fig.  18  Mass  flow  rate  in  regular  polygonal  ducts  with 
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Fig.  19  Shear  stress  distribution  along  outer  boundary  in  square  ducts 
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Fig.  22  Velocity  profiles  along  <t>  =  0  in  square  ducts  with  a  central 
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Fig.  23  Velocity  profiles  along  4>  =  tt/4  in  square  ducts  with  a  central 
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Fig.  24  Temperature  profiles  along  4>  =  0  in  square  ducts  with  a  central 
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Fig.  25  Temperature  profiles  along  4>  =  tt/4  in  square  ducts  with  a  central 
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Fig.  28  Temperature  profiles  along  <t>  =  0  in  20-sided  regular  polygonal  ducts 


Lra*  5'  j,sujfo€i.9^ni.€  bto^rj^a  sjoi;  0  «  0  Jv  SO**8^11?  TS'-  bojXdoi  1  *. 


ee 


<¥l  SYS) 


X  TO 


100 


£- 


01  x 


/>ra  ;'VsVqs'\ 


T  A 


101 


1 1  .-P)\(3'PJ 


•  •  • 

o  o 

n\(-2L  { -  Sl\ 


102 


Fig. 31  Nusselt  number  results  for  regular  polygonal  ducts 


with  a  central  circular  hole. 
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Fig.  32  Coordinate  system  for  circular  duct 


with  a  central  regular  polygonal  hole. 
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Fig.  33  Mass  flow  rate  in  circular  ducts  with  a  central  regular 
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Fig.  34  Shear  stress  distribution  along  inner  boundary  in  circular  ducts 
with  a  central  regular  polygonal  hole  for  7=2. 
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F ig o  37  Velocity  profiles  along  <t>  =  0  in  circular  ducts  with  a  central 
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Fig.  45  Coordinate  system  for  N-sided  regular  polygonal 


duct . 
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Fig.  46  Average  Nusselt  numbers  for  h yd r od ynami c a 1 1 y  and 
thermally  developed  laminar  flow  in  regular  poly¬ 


gonal  ducts. 
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Fig,  47  Temperature  profiles  along  =  tt/N  in  regular  polygonal  ducts  with 
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